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Abstract 

The impact of the non-2-spinon excitations of the one-dimensional S = 1/2 
Heisenberg antiferromagnet on the integrated intensity, the susceptibility, the 
frequency moments, and the Euclidian time representation of the T = dy- 
namic spin structure factor S{q, oj) is studied on the basis of finite-size data 

for chains with up to = 28 sites. 
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I. INTRODUCTION 



The spectral weight in the dynamic structure factor 

+00 

^''^ —00 

of the one- dimensional (ID) 5 = 1/2 Heisenberg antiferromagnet 

N 

H = jJ2Si-Si+i (2) 

1=1 

at T = has long been known to come predominantly from the 2-spinon excitations,0 which 
form a continuum between the lower and upper boundaries (with energies in units of J)i'0 

^L{q) = ^smq, u;[/(g) = 7rsin|. (3) 

However, the exact 2-spinon part of S{q,uj), here named S^'^\q,u!), was determined only 
very recently,i'@ based on new advances in quantum groups.i 

The exact expression of S^'^\q,uj), which had a fairly complex structure in its original 
form,i'i was simplified considerably in Ref. ^, evaluated numerically, and plotted versus q 
and uj. The spectral-weight distribution was found to diverge at ujL{q), 



S^'\q,^) - ^=^—J\n (4) 

and to vanish in a square-root cusp at uju{,q), 

S^'^\q,u) ~ \JuJu{q) - UJ. (5) 
These singularities differ from those of the approximate expression 

for the 2-spinon dynamic structure factor, which had been inferred from finite- data, sum 
rules, and Bethe-ansatz calculations.0 Otherwise, the line shapes of S^^\q,uj) and S^"'\q,uj) 
do not differ all that much.il 



The 2-spinon excitations were shown to account for 82.18% of the total intensity in 
S{q, uj). Where in (g, co')-space is the remaining spectral weight, from which classes of Bethe- 
ansatz solutions does it originate, and how does it affect various quantities that can be 
derived from the dynamic structure factor? 

A study of the 4-spinon contribution to S{q, u) along the lines of Refs. will shed 
light on the first two questions.!^ In answer to the third question, we investigate here the 
effects of the non-2-spinon excitations on four quantities which are related to S{q, u) and 
which can be computed with high precision from finite-A^ data for the ground-state wave 
function. 



II. INTEGRATED INTENSITY 

The integrated intensity of S{q,u!), i.e. the static spin structure factor 

I^Siq,co), (7) 


has been determined with high precision for wave numbers q < ISn/lA from finite- A^ data 
of cyclic chains with A^ < 28 sites.i This result is plotted in Fig. 1 for comparison with 
the exact 2-spinon integrated intensity I^'^\q) calculated from S^'^\q,to) via (0) and the 
integrated intensity0 

&\q) = ^ln^-±^^ (8) 
2n cos(g/2) ^ ^ 

obtained from the approximate result S^°'\q,uj). 

The finite- A^ data indicate that /(g) increases linearly from zero for small q and diverges 
logarithmically at g = vr. The initial rise of the finite- A^ data, /(g) 0.271g/7r,i is signif- 
icantly steeper than that of the 2-spinon contribution, I^'^\q) 0.237g/7r and that of the 
approximate result, I^°'\q) — > 0.25g/7r. Hence the integrated intensity of the non- 2-spinon 
part of S{q,u!) increases linearly in g too. 

The function /'■'^^(g) approximates the 2-spinon integrated intensity I'^'^^q) quite well for 
g/TT ^ 0.6. At larger g, this is no longer the case. The divergence predicted by expression 



(§), I^""^ ~ — ln(l — g/vr), is weaker than the divergence of the exact 2-spinon result, I^'^\q) ~ 
[— ln(l — Q'/tt)]'^/^. The inset of Fig. 1 shows the relative non-2-spinon integrated intensity, 
A/(2)(g) = 1 - /(2)(g)//(g), and the relative deviation, ^&\q) = 1 - I^^^q)/ I{q), of the 
approximate result If it can be assumed that the leading singularity of I{q) at g = tt is 
produced entirely by the 2-spinon part of S{q,uj), then the function AI^'^\q) must approach 
zero as g — > TT. The dashed line in the inset does not rule out that this assumption is correct. 
Whether or not the 4-spinon excitations contribute to the leading singularity in /(g) remains 
to be seen.i 

It is interesting to compare these results with the exact integrated intensity of the 
Haldane-Shastry model,i I^^^\q) = — (l/4)ln(l — g/vr), where non-2-spinon excitations 
have zero spectral weight in S{q,uj). It turns out that for g < 137r/14, I^^^\q) is a much 
better approximation of /(g) than I^'^\q) is.0 



III. SUSCEPTIBILITY 



The g-dependent susceptibility at T = is related to S'(g, a;) via the sum rule0 

oo 

= - /— %,^). (9) 



This quantity, which has been determined with considerable accuracy from finite- data, 
is plotted in Fig. 2 for comparison with the exact 2-spinon susceptibility x*^^"* {l) calculatedB 
from S^'^\q,ijj) via (|^) and the approximate resultS 

x^-^g) = (10) 

vr"^ sm g 

inferred from (|^). 

The normalization of S^'^\q^uj) was chosen such that the exact value of the direct 
susceptibilty,0 x(0) = I/tt^, is correctly reproduced. With increasing g, x^"'\l) deviates 
in a downward direction from x{q)- Its divergence at g = tt, x^''\q) ~ (tt — g)^^, is slightly 
weaker than that of the exact 2-spinon susceptibilty,@ ~ \/~ ~ q) / which, 

presumably coincides with the leading singularity of xio)- 



The contribution of the non-2-spinon spectral weight of S{q, u) to in the hmit g ^ 
is small but non-negligible as indicated by the result 7r^x'^^^(0) = 0.9500..., which falls some 
5% short of the exact direct susceptibilty. The relative non-2-spinon contribution to the 
susceptibility, Ax*^^)(g) = 1 — X^'^Kq) / xi^) ^ stays smaller for all g- values than the relative 
non-2-spinon integrated intensity AI^'^\q) (see inset). This indicates that the non-2-spinon 
spectral weight is located predominatly above the 2-spinon continuum. 



IV. FREQUENCY MOMENTS 

Yet a different way to assess the non- 2-spinon part of S{q,u) employs the frequency 
moments 

K4q)^ f ^u^S{q,u), n = l,2,..., (11) 



which are related, via sum rules,0 to short-range multi-spin correlations in the ground state. 
For n = 1 we know the exact results,^ 

i^i(g) = |§(l-cosg), (12) 

where Eq = — A^(ln2 — 1/4) is the ground-state energy. For n = 2,3,4,5 high-precision 
results have been calculated from finite- iV data for the associated ground-state expectation 
values.! The moments K^\q) of the exact 2-spinon dynamic structure factor S^^''{q,u) have 
been determined in Ref. § and the moments K^f'^q) of S'^"'\q,uj) in Ref. |12|. 

For n = 1, both Ki'\q) and Ki'\q) reproduce the g-dependence of the exact sum rule 
(|l^) correctly, but the prefactors are smaller, 

= 0.8039..., = 0.8462..., 13 

which again reflects the missing spectral weight of the non-2-spinon excitations. The q- 
dependence of the moment ratios 

i?,(g)^:^, forn = 2,3,4,5, (14) 
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of the full dynamic structure factor as inferred from finite- data are shown in Fig. 3 along 
with the corresponding moment ratios R^\q) of S^'^\q^uj) and the moment ratios R^\q) 
of S^'^\q,uj). The most striking observation is that Rn{q) approaches a nonzero value as 
g — > 0, whereas the exact and the approximate 2-spinon moment ratios both go to zero: 
^uHq) ~ ^I^Kq) ~ q^~^- This means that for long wavelengths the frequency moments 
Kn{q)-,n > 2, are dominated by non-2-spinon excitations, which are necessarily located 
above the narrow 2-spinon band. In other words, the 2-spinon dynamic structure factor 
S^'^\q,uj) does not contribute to the leading O(g^) term of Kn{q) for n>2. 

At larger wave numbers, the impact of the non- 2-spinon excitations on the moment ratios 
is more modest but still significant. Here the deviation of R^\q) from Rn{q) is almost q- 
independent, and it grows with increasing n. This again indicates that the non- 2-spinon 
spectral weight comes for the most part from higher frequencies than the 2-spinon spectral 
weight. 

The moment ratios R^^^q) agree very well with R^\q) at small g, but then deviate 
upwardly. For q ^ 7r/3, they rise even above the ratios Rn{q)- This discrepancy, which 
becomes more conspicuous with increasing n, reflects the fact that S^'^\q, u) underestimates 
the spectral weight near the lower continuum boundary UL^q) and overestimates the spectral 
weight near the upper boundary uju{q). 



In order to study the significance of the non- 2-spinon excitations in a dynamical quantity 
by the same kind of comparison, it is useful to consider the Laplace transform of the dynamic 
structure factor,0 





which can be interpreted as a Euclidian time representation of S{q, uj). For r = 0, this is the 
integrated intensity (0). From finite- iV data for S{q,uj) as obtained via recursion method 



V. EUCLIDIAN TIME REPRESENTATION 




(15) 
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for systems with N <28 sites,cJ this quantity can be accurately extrapolated to ^ oo if 
q 71. For the graphical representation, it is convenient to plot the function 

p{<l,t) = ^Jf'^l - 1, t = Juji{p)Texp[uJi{p)T], (16) 
S{q,T) 

which was used for the finite- iV extrapolation, instead of S{q, r) itself. The resulting curves, 
shown as solid lines in Fig. 4 for several g-values, rise from zero with zero initial slope and 
then become almost linear in t with a g-dependent slope. 

If the threshold singularity were a square-root divergence as predicted by (P), then the 
asymptotic growth of p{q, t) would be exactly linear. The logarithmic correction in the exact 
2-spinon threshold singularity (^), however, leads to a slight modification of the asymptotic 
growth of p(g, t), 

The dashed lines in Fig. 4 show the function p^^^ (g, t) as inferred from the exact 2-spinon 
dynamic structure factor. The discrepancies are fairly small over the range of t shown. The 
deviation A^'^\q,t) = p^'^\q,t) — p{q,t) is shown in the inset. Not surprisingly, the function 
p^"'\q,t) inferred from (^) deviates more strongly from p(q,t). This comparison was already 
made in Ref. 
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FIGURES 

FIG. 1. Integrated intensities (|^ in comparison: I[q) is represented by finite- data for 
iV = 6,8,...,28. is the exact 2-spinon result. I^"'\q) is the approximate result (|8|). The 

inset shows AI^^\q) = 1 - I^'^\q)/I{q) and Al^'^^q) = 1 - I^^^Xq) / 1 {q) . 

FIG. 2. The (/-dependent susceptibilty (^) in comparison: xio) is represented by finite-A^ data 
for N = 6, 8,. ..,28. x^'^Ko) is the exact 2-spinon result. x^^H?) is the approximate result (p!oD. 
The inset shows the relative non- 2-spinon integrated intensity Ax*-^^(g) = 1 — X^'^\l)/x{l) 
Ax(")(g) = l-x("ng)/x(g). 



FIG. 3. Ratios of frequency moments (|14D. For n = 2,3,4,5, the Rnio) represent finite- A'^ 
data for N = 6, ...,28. The Rn\q) are exact 2-spinon results and the Rn\q) are the moment 
ratios for (P). 

FIG. 4. The function 1^ for q = tt/A, 7r/3, 7r/2, 27r/3, 37r/4. (bottom to top) The solid lines 
represent extrapolated finite- data and the dashed lines the exact 2-spinon-part of that function. 
The inset shows the non- 2-spinon part A(2)(g,t) of the same function. 
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